Kerov's interlacing sequences and random matrices 
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q Abstract 

To a N x N real symmetric matrix Kerov assigns a piecewise lin- 

j» ear function whose local minima are the eigenvalues of this matrix and 

whose local maxima are the eigenvalues of its (N— 1) x (N — l) subma- 
trix. We study the scaling limit of Kerov's piecewise linear functions 
for Wigner and Wishart matrices. For Wigner matrices the scaling 
limit is given by the Verhik-Kerov-Logan-Shepp curve which is known 
from asymptotic representation theory. For Wishart matrices the scal- 
ing limit is also explicitly found, and we explain its relation to the 
Marchenko-Pastur limit spectral law. 

1 1 Introduction 

Consider two sequences of real numbers {xi}f =1 , {yj}™z\ such that 

in 

^ x\ > yi > x 2 > ■ ■ ■ > x„-i > j/n-i > x n . (1.1) 

We say that the sequences {xi}f =1 , {yj}" =1 interlace. Define 



PL, 



(N 



n n—l 
i=l j=l 



X 

Following Kerov (see [5]) we define a rectangular Young diagram w^ Xi ^'^ v ^ (x) 
which is uniquely determined by the following conditions (see Figure [TJ : 

1) lo^Mty}^) : ]R — ). M is a continuous piecewise linear function, and 

d 

w { x i},{yj} ( x ) = -M except for finitely many points, which are exactly the 

ox 

local extrema of the function u/i^bw/^), 

2) {xi}f =1 are local minima of w^ Xi ^^ y ^ (x), {yj^Zi are local maxima 
of w* Xi }'* v >}(x), and there are no other local extrema. 

3) vj{ Xi *'{ y i}(x) = \x — Zq\ when \x\ is large enough. 

Let S be a N x N real symmetric matrix. By S we denote its (N — 
1) x (N — 1) submatrix; it is obtained from S by removing the A^th row 
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Figure 1: A part of a rectangular Young diagram 



and column. It is well-known that the eigenvalues of S and S interlace (see, 
e.g, p. 185]). Thus, to any symmetric matrix we can assign a rectangular 
Young diagram built from the eigenvalues of S and S. 

Let {Zij}°°j =1 be a family of independent and identically distributed 
zero-mean, real-valued random variables such that EZj^ = 1 and 

E|Zn| fe <oo, for all k = 1,2,3,.... 

A symmetric N x N matrix Xn with entries 

XN(hj) = X N (j, i) = Zij, for i <j, 

is called a Wigner matrix. Let w£(x) be the rectangular Young diagram 
constructed from the eigenvalues of and X^. Note that wj^(x) is a 
random function. We are interested in the limit behaviour of w£(x) as 
N oo. 
Let 

xarcsin(|) + \/4 — x 2 J , \x\ < 2, 

\x\ > 2, 

be the Vershik-Kerov-Logan-Shepp curve (see [T2] and [TU]). 
Theorem 1. As N — > oo ; we have 
1 



n(x) 




lim sup 



w§(xVN) -n( x ) 



0, 



in probability. 
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Figure 2: The rectangular Young diagram wn{x) and fl(x) for iV = 1100 



Interlacing sequences arise naturally in several areas. They provide a 
useful parametrization of Young diagrams (see [9] , [1] ) . They also appear as 
roots of two consecutive orthogonal polynomials (see [6j). A more general 
notion of interlacing measures was studied in [S]. 

For the first time the curve Q(x) appeared from a representation theo- 
retic problem. This curve is the limit shape of the random Young diagram 
distributed according to the Plancherel measure (see [12], [ID], [IB], and [H 
Section 5] for more details). Then it was found that tt(x) is a scaling limit 
for separation of roots of orthogonal polynomials (see [6]). This curve also 
arises as the scaling limit in the evolution of continuous Young diagrams 
(see [7]) and in random matrix theory. Let us formulate Kerov's result from 
[6J related to random matrix theory. 

Let h]\f d ~R N be a random hyperplane such that G hjy and the normal 
vector to Jin is uniformly distributed on the unit sphere. Let p be the 
projection operator to h^. Regard X^ as an operator in M. N and consider 
the operator pXjyp in hjy. The eigenvalues of Xn and pX^p interlace. Let 
us construct the rectangular Young diagram wn from these eigenvalues as 
above. 

Theorem ([BJ, Th. 3.6). As N — > oo, we have 



and the limit is uniform on x € M. 

Remark 1. In the context of Theorem Q] we consider the restriction to a fixed 
hyperplane while here the hyperplane is random. Another difference is that 




(1.2) 
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Theorem [I] establishes the convergence in probability while (1.2) gives only 
the convergence of mean. 

Let us proceed to Wishart matrices. Let M = M{N) be a sequence of 
positive integers such that 

M(N) 
lim — — = a > 1. 

N^oo N 

Let Wtv beaJVx M(N) matrix with i.i.d. entries of mean zero and variance 
1, and such that 

E\W N (l, l)\ k < oo, for all fc = 1,2,3, 

By Yn we denote the N x N Wishart matrix Yjq = WatW^. Let wjf(x) be 
the rectangular Young diagram which is defined by the interlacing eigenval- 
ues of Y/v and YXr. 

Let us define a continuous function Q : M — > R. Set 

= ^^f"^ = , x G [(a + 1) - 2^a; (a + 1) + 2^], 

and 



^a(^) = \x — a|, x S (— oo; (a + 1) — 2\/a] U [(a + 1) + 2\fa\ oo). 
These formulas determine the function £l a (x) uniquely. 
Theorem 2. As N — > oo, we have 



lim sup 



wJf(xN) — £l a (x) 



N 



0. 



in probability. 



Remark 2. The limit shapes £l a (x) are closely connected with Biane's limit 
shapes (see [2]). This connection is described in Section 5. 

This paper is organized as follows. In Section 2 we give some prelimi- 
naries. In Section 3 we prove Theorem [T] In Section 4 we prove Theorem 
[2] In Section 5 we describe a link between these limit shapes and semicircle 
and Marchenko-Pastur limit laws. 
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Figure 3: The rectangular Young diagram wj^(x) and Q a {x) for a = 2.25 
and N = 400 



2 Continual Young diagrams 

A continual Young diagram (see [5]) is a function w(x) on R such that 

1) |u>(a;i) — u>(x2)| < |^l — ^2! for any si, £2 £ R- 

2) There exists a point xq 6 R, called the center of ui, such that = 
|x — xo| when |x| is large enough. 

The set of all continual Young diagrams is denoted by T>. For any w E T> 
we define a function cr(x) 



Since cr(x) satisfies the Lipschitz condition 1), its derivative a'(x) exists al- 
most everywhere and satisfies |cr'(a;)| < 1. Note that cr'(x) is compactly 
supported. The function w(x) is uniquely determined by cr'(x). Moreover, 
w(x) is uniquely determined by the second derivative a"(x), which is under- 
stood in the sense of distribution theory. 

Let us define the function pk : V — > E, k S N, by setting 



a{x) = ~(w(x) 



x\). 




It is easy to see that for the rectangular Young diagram w 
(see Introduction) we have 
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We will need the following fact. 

Lemma 2.1 Q2], Lemma 5.7). Let J-([a;b]) be the set of all real-valued func- 
tions f(x) supported on the interval [a,b] G M. and satisfying the Lipschitz 
condition |/(a?i) — f(x2)\ < \xi — X2\- 

On the set J-([a;b]), the weak topology defined by the functionals 

/(s)-> / f{x)x k dx, fc = 0,1,2,..., 

Jxe[a;b] 

coincides with the uniform topology. 

3 Proof of Theorem [T] 

Lemma 3.1. Let be the (random) rectangular Young diagram defined 
by the eigenvalues of Wigner matrices Xn and Xn (see Introduction). Then 

Pk(w%) ^ jo, k is odd, 

Nk/2 N -*°° l (fc/2)l(fc/2)l » k is even, 
where the convergence is in probability. 

Proof. The proof of this Lemma closely follows the proof of Wigner's the- 
orem (see, e.g., [H Section 2.1]) and is based on the well-known moment 
method. 

Let {\f}fL 1 be the eigenvalues of X^ and let {A^" 1 }^ 1 be the eigen- 
values of X/v- We have 

. / N N-l \ 

p k (w§) = £ * (Af ) - £ 5 (Af - 1 ) \dx = tr (x%) -tr ( X %_ x 

Jxm y i=1 j=1 J 

Let i^v = (ii>*2i • ■ • >ifc) range over the set of multi-indices such that 1 < 
ii, i 2 , ■ ■ ■ , ik < N. Likewise, let ijv-i = (11,12, ■ ■ ■ ,ik) range over the set of 
multi-indices such that 1 < i\, i 2 , ■ . . , ik < N — 1. Then 

tr (X^j - tr (X N _^ = ^2x N {i 1 ,i 2 )X N (i 2 ,h) ...X(i k ,H) 

ijv 

- ^ X N (i 1 ,i 2 )X N (i 2 ,i 3 ) ■ ..X(i k ,ii) = -XjvCii^-Xjvfa,^) . . .X(i fc ,«i), 

ijv-i ijv:JVeijv 

where the last sum is taken over indices i^v = • • • > ik) such that there 

exists r, 1 < r < k, satisfying i r = N. 
First, we compute 

E ^2 X N (ii,i 2 )X N (i 2 ,i 3 ) . . . X{i k ,i\) • 

\iiv:Afeijv / 
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This sum can be written as a sum of terms corresponding to suitably defined 
graphs that are in their turn associated to words. Suppose k is odd; then 
the same estimates as in [U Lemma 2.1.6] show that the contribution to the 
degree n k is equal to 0. 

Suppose k is even; then the main contribution is given by the so called 
Wigner words (see [U Def. 2.1.10]). The number of Wigner words is equal to 
A;! 

— — r— — — -. — — . The only difference of our case with the case of Wigner's 

(k/2 + l)!(fc/2)! 

theorem is that one vertex of a graph should be labeled by the special symbol 
N. This condition gives an extra factor (k/2 + 1). Thus, we obtain 



lim 7V" fc / 2 E V X N {i 1 ,i 2 )X N (i 2 ,i 3 )...X(i k ,i l ) 

N^oo \ * — ' 

ijv^Sijv 



0, k is odd, 

fej . ■ (3-2) 

- — — — — — -, k is even. 
(*/2)!(fc/2)!' 



Secondly, we have 



lim \N- k / 2 V X N (i 1 ,i2)X N {i2,i 3 )...X(i k ,i 1 )) =0. (3.3) 



Indeed, this equality can be proved in the same way as in [TJ Lemma 2.1.7]. 



From (3.2) and (3.3) it follows that p k converges to the right-hand side 
of (3.1) in L 2 and, consequently, in probability. 

□ 



Lemma 3.2 Q4J Prop. 5.3). We have 
p fc (0) = < 



0, k is odd, 

k is even. 



l(k/2)l(k/2) 



Lemma 3.3. There exists an interval [-B; B] such that the probability that 
the inequalities 

-B < \% < ■ ■ ■ < Af < B 
hold tends to 1 as N —> oo. 

Proof. This is a well-known fact from random matrix theory, see, e.g. [Tj. 

□ 

Let z/v be the center of the rectangular Young diagram w^. It is clear 
that 

lim —p= = lim — - ' - =0, in probability. 
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This equality and Lemma 3.3 imply the convergence of to Q(x) outside 
the interval [— B;B]. 



(3.4) 



From Lemmas 3.1 and |3.2| for any k G N we have 



lim . , 

AM-oo N k ' 2 



pk(£l), in probability. 



The convergence of wjy to Q(x) inside the interval [-B; B] follows from (3.4) 
and Lemma |2. 11 



4 Proof of Theorem [2] 

Lemma 4.1. Let w N be the (random) rectangular Young diagram defined 
by the eigenvalues of Wishart matrices Y/v an d ^JV (see Introduction). Then 
the following limit exists 

— — > rrifc, in probability. (4-1) 

N K TV— >oo 

The generating function of rrik is given by the equation 

^ , \ -, k 1 I -, (a — l)z + l \ 

G a (z) :=1 + Y,m k z k = - 1 + 



k=l 



2 V y/(a-l) 2 z 2 -2(a + l)z + l 



Proof. The proof of this Lemma is based on moment method and follows 
that given in [H Exercise 2.1.18]. 

Let 'in = i%, . . . ,ijv), Jn = (jl,j2, ■ ■ • >ij\r) range over the set of multi- 
indices such that 1 < ii, Z2? • • • j*jv < 1 < ji>j2, ■ ■ ■ ,3n < M(N). We 
have 

Pfc(«&) = tr(Y#) - tr(y^_j) = Kjv-(«i J *2)i r Jv-(*2,*3).--yiv(»j>r,n) 

ijv:Afeijv 

= ^ VWjv(«i, ji)VWjv(i2, ii)Wiv(«2, ji)yVjv(?2, j 2 ) • • • Wjvfe,ifc)Wjv(n, jit), 

ijvjjv^eijv 

where the condition N £ means, as before, that there exists r such that 
i r = N. 

The limit relation 

lim — J! = 

N^oo N 2k 

can be proved in the same way as in [TJ Section 2.1]. Thus, the only thing 
we need is to compute the leading term in 

E W N {ii,j 1 )W N {i 2 ,ji)W N (i2,ji)W N (i2,j2)---WN(ik,jk)W N {i 1 ,j k ) 

yiwjjv^eijv 

(4.2) 
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Let us recall that a Dyck path T>\ of length 21 is an integer-valued se- 
quence {S n } < n <2i such that S = 0, S 2 i = 0, \S(i) — S(i — 1)| = 1 for 
1 < i < 21, and > for < i < 21. Let a(T>i) be the number of indices 
i such that S(i) — S(i — 1) = — 1 and S(i — 1) is odd and let b(V{) be the 
number of indices i such that S{£) — S(i — 1) = — 1 and S(i — 1) is even. We 
have a(X>j) + = Z. 

There is a bijective correspondence between Wigner words and Dyck 
paths (see [U Section 2.1]). It is easy to see that the main contribution to 
(4.2 ) has order N k and is given by Wigner words of length 2k or, equivalently, 
by Dyck paths of length 2k. If the sequence 12, 32, ■ ■ ■ , ik, jk) is a 

Wigner word then a(T> k ) is equal to the number of distinct j-indices in the 
tuple and b(T>k) + 1 is equal to the number of distinct i-indices in the tuple, 
where T> k is the corresponding Dyck path. Therefore, each word gives a 
contribution (b(V k ) + l)a a(x,fe) and we obtain 

m k = (b(V k ) + l)a a ^\ (4.3) 
all v k 

Let ft be a formal variable, and let 

d r ■= P r ~ a[Vr) a a{Vr \ r > 1, d = 1; 
all v r 

e r ■= ft a ^a r ~ a{Vr) , r > 1, e = 1. 
all x> r 

Considering the moment of the first return of Dyck path D r to zero we 
obtain 

r r 

dy Gi dy — j Gj — \ , (3y> ft — j^j — 1 ' 

3=1 3=1 

Hence, 

a 

We obtain 

r 

d r = ft ^2 d r -jdj-i + ad r -\. (4-4) 
3=2 

Let d(z) be the generating function of {d r }: 

oo 

r 



:= 1 + d 



Using (4.4), we get 

d( z ) = 1 + /3zd(z) 2 + (a - ft)zd(z). 
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Solving this equation and choosing the sign from the condition d(Q) = 1, we 
obtain 

l-(a-P)z-y/((a-P)z-l)2-4l3z 
d{2) = Wz • 



From (4.3 ), we get 



d 



k=l P 



This completes the proof of Lemma 4.1 



□ 



For a > 1 let fi a (cc) be defined by the formula 
1 



f2 a (x) '■= < 



7T 



(2a — x) arcsin 



q+l— x 
2J^c 



arctan 



{a-l) 2 ~x(a+l) 



+ yj2a + 2x-l- (x-a) 2 



(a-l)V \x-a) 2 +2a+2x-l 

for x G [a + 1 - 2y/a; a + 1 + 2 



Jx — a|, for x G (— oo; (a + 1) — 2-y/a] U [(a + 1) + 2yfa; oo), 
and for a = 1 let 



fti(x) := 



f i f (x - 2) arcsin(^x - 1) - ^4 - (s - 2) 2 ) + |, 
(jx-l|, for x G (— oo; 0] U [2; oo). 



for x G [0; 2] 



It is readily seen that for a > 1, f2 Q (x) is a continual Young diagram 
with the center a. 

Lemma 4.2. The generating function of {pk(^a)} is given by the formula 



k=l 

Proof. We recall that 



1 + 



(a-l)z + l 



2 I ^/(a-l) 2 z 2 -2(a + l);z + l 7 ' 



dx, k G N. 



It can be verified by a direct computation that for a > 1 

x G [q + 1 - 2^/a; a + 1 + 2^/a], 



TO = 



x + (a — 1) 
■KXyfAa. — (x — (a + l)) 2 ' 
0, x^[a+l-2v / a;a + l + 2 v /a]. 



Using the technique of Stieltjes transform (see, e.g., [U Section 2.4]) we 
obtain the statement of the Lemma. □ 
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Let zjj be the center of the diagram wjj. Note that 
4 Y N (N,N) Z?=iW(N,j)W(j,N) . 

-w = n = n ^ 111 p robabiht y- 

This fact implies that the center of jrwJj(Nx) converges to the center of 



Prom Lemmas 4.1 and 4.2 we have 



— ^ir ► Pki^a), m probability. 

Also it is known that the spectrum of Wishart matrix is supported by a 
fixed compact interval with probability close to 1. Theorem [2] follows from 
these facts in the same way as in the proof of Theorem [T] 

5 Connection with semicircle and Marchenko-Pastur 
laws 

In this Section we briefly describe a connection between limiting shapes O, 
fl a and well-known semicircle and Marchenko-Pastur distributions, respec- 
tively. 

For an interval / by we denote the set of probability measures 

which are supported by /. For a measure fj, £ M.(I) let 

:= J x k dfi, k G N. 

By T>{I) we denote the set of continual Young diagrams such that cr"(x) is 
supported by /. 

Lemma 5.1. There is a bijective correspondence \x — > w between A4(I) and 
It is characterized by the relation 

1 + f>z fc = exp ^ A (5.1) 

fc=l \k=l / 

Proof. See [5], [8] . □ 

The measure \i is called the transition measure of w. For more details 
about this correspondence, see [5], [H Section 8]. 



A direct computation of the left- and right-hand sides of (5.1) leads to 
the following Proposition. 

Proposition 5.1. a) The semicircle law is the transition measure ofQ. 

b) The Marchenko-Pastur distribution with parameter a is the transition 
measure of Vt a . 
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The item a) was first noted in [5]. 

In [2] Biane introduced a family of curves which appeared as scaling lim- 
its in some problem of asymptotic representation theory. As shown in 
the transition measure of every Biane's curve coincide with a Marchenko- 
Pastur distribution, within a homothetic transform. Thus the Biane's curves 
and the curves Q a are closely related. 

Remark 3. Let x n = {x\, X2, ■ ■ ■ , x n ) be a sequence of real numbers and 
suppose that for every n the sequences x n and x n _i interlace. Let /i be 
a probability measure on R and let 5{x) be the Dirac measure at i 6 I. 
Assume that 



It was shown in [6] that there are sequences with different limiting measures 
fi but with the same scaling limit of y; x "' x «-i. Therefore, the convergence of 
n Y^i=i $( x i) to the measure \i does not imply the convergence of u; x ™> x ".-i 
to the continual Young diagram with the transition measure fi. 
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